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Abstract 

Chemoautotrophy has been little studied in typical coastal marine sediments, but may be an important component of 
carbon recycling as intense anaerobic mineralization processes in these sediments lead to accumulation of high amounts of 
reduced compounds, such as sulfides and ammonium. We studied chemoautotrophy by measuring dark-fixation of ^^C- 
bicarbonate into phospholipid derived fatty acid (PLFA) biomarkers at two coastal sediment sites with contrasting sulfur 
chemistry in the Eastern Scheldt estuary, the Netherlands. At one site where free sulfide accumulated in the pore water right 
to the top of the sediment, PLFA labeling was restricted to compounds typically found in sulfur and ammonium oxidizing 
bacteria. At the other site, with no detectable free sulfide in the pore water, a very different PLFA labeling pattern was found 
with high amounts of label in branched i- and a-PLFA besides the typical compounds for sulfur and ammonium oxidizing 
bacteria. This suggests that other types of chemoautotrophic bacteria were also active, most likely Deltaproteobacteria 
related to sulfate reducers. Maximum rates of chemoautotrophy were detected in first 1 to 2 centimeters of both sediments 
and chemosynthetic biomass production was high ranging from 3 to 36 mmol C m~^ d~\ Average dark carbon fixation to 
sediment oxygen uptake ratios were 0.22±0.07 mol C (mol 02)~\ which is in the range of the maximum growth yields 
reported for sulfur oxidizing bacteria indicating highly efficient growth. Chemoautotrophic biomass production was similar 
to carbon mineralization rates in the top of the free sulfide site, suggesting that chemoautotrophic bacteria could play a 
crucial role in the microbial food web and labeling in eukaryotic poly-unsaturated PLFA was indeed detectable. Our study 
shows that dark carbon fixation by chemoautotrophic bacteria is a major process in the carbon cycle of coastal sediments, 
and should therefore receive more attention in future studies on sediment biogeochemistry and microbial ecology. 
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Introduction 

Reoxidation of reduced intermediates like sulfide and ammo- 
nium formed during anaerobic mineralization processes is an 
important process in coastal marine sediments. Oxygen is typically 
only found in the top millimeters of these sediments and along 
macrofauna burrows [1], and carbon mineralization proceeds in 
general by anaerobic processes primarily sulfate reduction. This 
results in the production and accumulation of large amounts of 
reduced compounds such as various forms of reduced sulfur and 
ammonium [2]. In typical coastal sediments, free sulfide in the 
porewater is however often only detected below a couple of 
centimeters as it quickly reacts with iron hydroxides forming iron 
sulfide (FeS) or pyrite (FeS2) [3] . Only in very active sediments or 
sediments containing little reactive iron, free sulfide can be found 
near the oxic top layer [3]. Long term burial of reduced 
compounds is thought to be a minor process [3] and they are 
mostly transported to more oxidized horizons by either diffusion 
or bioturbation [4] . Oxygen is eventually the main oxidant of these 
reduced compounds although intermediate reoxidation steps by a 
variety of anaerobic pathways using nitrate or iron and manganese 



oxides may also be important [3]. It is estimated that reoxidation 
processes on average explain 70% of the sediment oxygen flux in 
shelf sediments [5] and this value is expected to be higher in active 
intertidal areas as anaerobic mineralization will be more 
important. 

Many of the known prokaryotes involved in reoxidation 
processes are chemo(litho)autotrophs that use the energy gained 
from inorganic reactions to grow by fixing inorganic carbon in the 
dark [6] . Chemoautotrophic carbon fixation has been shown to be 
an important process in, for instance, extreme marine ecosystems 
such as hydrothermal vents [7,8] and in the chemocline of anoxic 
marine basins [9,10]. The current consensus is however that 
chemoautotrophy is a relatively minor process in coastal sediments 
due to the relatively low growth yields of chemoautotrophic 
organisms and the competition with chemical oxidation reactions 
[3] . In addition, true chemoautotrophic bacteria have to compete 
with mixotrophic and heterotrophic bacteria that are able to 
oxidize reduced sulfur compounds [11], which could be relevant 
especially in active coastal sediments receiving large amounts of 
organic matter. Studies where chemoautotrophy was actually 
quantified by determining dark carbon frxation rates are rare for 
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typical coastal marine sediments and we have only been able to 
locate four studies: two on shallow subtidal sediments from the 
Baltic [12,13], one study on an intertidal sand flat in the German 
Wadden Sea [14] and a recent study on three brackish coastal lake 
sediments in Brazil [15]. However, recent estimates suggest that up 
to 0.29 Pg C y could be potentially fixed by chemoautotrophic 
microorganisms in near shore and shelf sediments worldwide 
compared to 0.92 Pg C y~' of mineralization [16], suggesting a 
major role in the sediment carbon cycle. Finally, the dominant 
chemoautotrophic bacteria involved in sulfur oxidation are not 
well known in coastal marine sediments. A recent study identified 
an uncultered group of Gammaproteobacteria as important 
players [14], but there may be many other groups involved in 
the diversity of reoxidation processes that occur in marine 
sediments. 

We studied chemoautotrophy in two intertidal sites with 
contrasting sulfur chemistry: a site where free sulfide was not 
detected in the top few centimeters of the sediment and a very 
active site where high concentrations of free sulfide were found 
right to tlie top of the sediment. The main substrates driving 
chemoautotrophy are therefore expected to be different at both 
sites, namely free sulfide at the very active site versus iron sulfides 
in the more typical coastal sediment. Chemoautotrophy rates were 
determined by incubating sediment cores with stable isotope 
labeled '^C-bicarbonate and measuring labeling in phospholipid 
derived fatty acids (PLFA). This method both yields estimates of 
total chemoautotrophy rates and provides an indication of the 
active bacterial community [17-19]. The diversity of Rubisco 
genes was studied to further indicate possible active chemoauto- 
trophs that use the Calvin cycle for carbon fixation. Finally 
chemoautotrophy rates were compared with diffusive oxygen 
fluxes and carbon mineralization rates. 

Materials and Methods 

Description of field sites 

Two field sites in the Eastern Scheldt estuary (The Netherlands) 
were selected, which were expected to show high mineralization 
rates and have major differences in sulfur chemistry. The site in 
the Zandkreek area (51°32'41"N, 3°53'22"E) was situated next to 
a Pacific oyster [Crassostrea ^as) bed and was sampled in April 
2005 (abbreviation ZK05) and October 2007 (ZK07). The Pacific 
oyster is an invasive species in the area that was introduced in the 
Eastern Scheldt around 1970. It stimulates sedimentation and 
sediment carbon mineralization either by decreasing water 
currents over the sediment or via pseudo-feces production and 
biodeposition [20]. Sediments were non-suUidic in the top 5 
centimeters in 2005 and slightly sulfidic below 2 centimeter in 
2007 (See Result). 

The Rattekaai site (51°26'2r'N, 4° 10' 11 "E) was situated at the 
entrance of a salt marsh creek where macroalgal debris (mainly 
Uba. derived) accumulates and is buried during winter. The 
sediment was highly sulfidic right to the top and samples were- 
taken from patches where the sediment was t:()vered with a whitish 
layer in April 2005 (RK05) and May 2006 (RK()6). Based on 
microscopy, typical Beggiatoa-like sulfur-oxidizing bacteria were 
abundant in the top few millimeters of the Rattekaai sediment, 
especially in 2005 and to a lesser degree in 2006. 

Ethics statement 

The Eastern Scheldt estuary is a Natura 2000 protected area 
and research permits for both sites were granted by the 
"Vereeniging Natuurmonumenten" and the Provence of Zeeland. 



Sediment sampling 

Undisturbed sediments were sampled with two sizes of 
polycarbonate core liners. The smaller cores (internal diameter 
4.6 cm) contained silicon-filled injection ports at every 0.5 
centimeter and were used for measuring chemoautotrophy rates. 
The larger cores (internal diameter 6 cm) were used for additional 
measurements of porewater profiles and sediment characteristic, 
and for measuring mineralization rates. Sediments were sampled 
at low tide and therefore did not have overlying water. Core's were 
processed the same day for chemoautotrophy rate measurements 
and other analyses. 

Chemoautotrophy rates 

Chemoautotrophy rate measurements were started by injecting 
100 Ml of 20 mM NaH'-^COa (99% ''^C; Cambridge Isotope 
Laboratories, Andover, MA, USA) horizontally into the sediment 
cores at 0.5 cm depth intervals by using the line-injection method 
[2]. The '^C-label was dissolved in artificial seawater lacking 
calcium or magnesium in order to avoid precipitation [21]. The 
label was made oxygen free by bubbling with nitrogen gas shortiy 
before injec tion. Sediment cores were incubated in the dark within 
2°C of the in-situ temperature (see Table 1) for various periods of 
up to 4 days, and were ventilated daily by removing the top 
stopper for one minute (ZK) or incubated without top stoppers 
(RK) to circum\'ent the development of suboxic condition in the 
headspaci-. Aftc-r incubation, sediment cores were sliced to a depth 
of 5 cm and sediment slices were quickly centrifuged (4500 rpm, 
5 min) to collect porewater for concentration and '^C analysis of 
dissolved inorganic carbon (DIC). Sediments were subsequenfly 
frozen at — 20 C and lyophHized before further analysis. 
Unlabelled, control cores were also processed. 

PLFA analysis and calculation of chemoautotrophy rates 

LyophUized sediments were analyzed for PLFA concentrations 
and C-labeling as described before [22,23]. In short, PLFA were 
extracted according to standard protocols and were analyzed by 
gas chromatography - isotope ratio mass spectrometry (GC- 
IRMS, Thermo, Bremen, Germany) on an a-polar analytical 
column (HP5-MS, Agilent, Santa Clara, CA, USA). Stable carbon 
isotope ratios are reported as 5'^C ratios on the VPDB scale. 
Excess C in individual PLFA was calculated as in Boschker et al 
[23] and divided by the atom percent excess '"^C in the DIC pool 
to calculate actual PLFA synthesis rates. Only very minor labeling 
was found in poly-unsaturated PLFA typical for Eukarya (see 
Results) suggesting that PLFA labeling was primarily by Bacteria. 
We therefore used the labeling data for all common bacterial 
PLFA in the 12:0 to 20:0 range in our calculations and not just the 
specific bacterial biomarker PLFA [24]. Total bacterial chemo- 
autotrophy rates were determined by summing synthesis rates in 
all PLFA typically found in bacteria and converted to chemoau- 
totrophic biomass production by dividing by the typical PLFA 
content of aerobic bacteria (55 mmol PLFA-C (mol biomass C)~' 
[24,25]). To study the diflerences in active chemoautrophic 
bacterial communities, we performed a principle component 
analysis (PCA) on log-transformed PLFA '^C-labeling data (in 
Mol%) using the Statistica software package (StatSoft, Tulsa, 
USA). 

Additional measurements 

Oxygen profiles were determined with oxygen microelectrodes 
(Unisense Ox 100, Aarhus, Denmark), which were lowered with a 
micromanipulator into the sediment until no oxygen was detected. 
Two profiles were recorded for each duplicate sediment core (four 
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Table 1. Sediment in-situ temperature, sediment characteristics, oxygen consumption rates, carbon mineralization rates, 
chemoautotrophy rates and yields (averages ± standard deviations, N = 2) for the coastal marine sediments in this study. 













O2 penetration 


O2 flux 


C mineralization^ 


Chemoautotrophy^ 


Yield C/O2 


Site/Year 


Temp. °C 


POC' (%) 




depth (mm) 


(mmol m^ d^; 


(mmol m^ d^) 


(mmol C n 


^^ d 


(mol C (mol Oj)^) 


RK05 


14 






0.45 ±0.10 


17.2±3.0 




5.5 ±1.9 




0.32 ±0.11 


RK06 


17 


2.0 


10.9 


0.23 ±0.06 


192±41 


197±36 


36.3 ±4.8 




0.1 9 ±0.03 


ZK05 


14 






1.7±0.1 


15.0±0.4 




2.6±0.3 




0.1 7± 0.02 


ZK07 


13 


0.6 


7.7 


0.95 ±0.06 


15.5±1.6 


105.9±19.1 


2.9±0.2 




0.18±0.01 



^Data for 0-1 cm sediment depth. ^ Data integrated over 0-5 cm sediment depth. 
doi:l 0.1 371 /journal.pone.Ol 01 443.t001 



profiles total), which were kept within 2°C of the in-situ 
temperature. Oxygen fluxes into the sediment were calculated as 
described in Van Frausum et al. [26] with sediment tortuosity 
estimated from sediment porosity as in Boudreau and Meysman 
[27]. 

Sediment porewater was sampled by shcing duplicate sediment 
cores in an anaerobic glove-box filled with 3% hydrogen in 
nitrogen gas (Coy Laboratory Products, Ann Arbor, MI, USA) 
and slices were centrifuged at 4500 rpm for 10 min at in situ 
temperature. Samples for sulfide analysis were immediately fixed 
in zinc acetate and analyzed according to Cline [28] . Samples for 
ammonium and anion analysis were frozen, and analyzed on a 
QuAAtro segmented flow analyzer (Seal Analytical, Norderstedt, 
Germany) and suppressed high performance ion chromatography 
on a Dionex lonpac AS- 14 column (Thermo, Sunnyvale, CA, 
USA), respectively. Samples for ' ^C-DIC were added to headspace 
vials (10 ml) and after acidification analyzed for DIG concentra- 
tions and ''^G-content by elemental analyzer - IRMS [29]. 

Sediment carbon mineralization rates were determined using 
the jar method [30] . Sediment cores were shced as above and were 
incubated in completely filled centrifuge tubes. Gentrifuge tubes 
containing the sediment were sealed in air-tight incubation bags 
filled with nitrogen gas to keep them strictly anaerobic and were 
incubated within 2°G of the in situ temperature for up to 6 days. 
Sediment porewater was collected and analyzed as described 
above. Mineralization rates were calculated from DIG and 
ammonium production with time and ammonium production 
was converted to carbon mineralization rates by using the 
sediment G/N ratio (Table 1). 

Rubisco type lA clone libraries 

To further study the diversity of chemoauto trophic bacteria that 
utilize the Galvin cycle for carbon dioxide fixation, Rubisco clone 
libraries were constructed for both sites in March 2008. Sediments 
were sampled as described above, and the top 0.5 cm of the cores 
showing the highest chemoautotrophy rates was collected and 
immediately frozen at — 80°G. Total community DNA was 
extracted from 0.5 g of wet sediment using the MoBio UltraGlean 
Soil DNA Isolation kit according to protocol (MoBio, Garlsbad, 
GA, USA). 

We developed a new degenerative primer set to specifically 
amplify Rubisco type lA as this group contains most of the true 
chemoautotrophic bacteria involved in sulfur and ammonium 
oxidation [31,32]. The new primer set also targets Begff,atoa hke 
Rubisco sequences [33], which was important as iJcggwtoa-like 
bateria were found at the RK site but were not covered by 
previously pubhshed primer sets developed for chemoautotrophic 
bacteria. The primer set also targets some of the lower branching 
Type IB sequences found in unicellular cyanobacteria, and 



consists of forward primer 571 (GAYTTYAGGAARGAYGAYG) 
and reversed primer 898E (AGRGGGAARTGRATRGG). The 
primer set was first tested against a positive control [Uiioalk- 
alimkrobium aerophilu kindly provided by Gerhard Muyzer, Delft 
Technical University, The Netherlands) and PGR conditions were 
subsequently optimized to specifically amplify the target sequences 
from sediment DNA extracts. 

The final PGR reaction mixture contained: 2.5 \il lOx standard 
Taq reaction buffer (without Mg), 3.0 mmol L~' Mg^"*^, 0.2 mmol 
L"' dNTPS, 0.2 ixmol L"' of each primer (571 and 898E), 2 U of 
NEB Taq DNA polymerase, 5% v/v DMSO, 0.2% w/v BSA and 
16 nl autoclaved demi water. The PGR cycling intervals were 
established as follows: preheating at 94°G for 5 minutes, followed 
by 40 cycles of denaturation step at 94°G for 1 minute, annealing 
step at 51°G for 30 seconds and extension at 72°G for 30 seconds. 
The PGR reaction was finished with a final extension time of 7 
minutes at 72°G. PGR products for each sample (RK and ZK) 
were cloned into Escherichia coli Top 10 cells using TOPO TA 
cloning kit (Invitrogen, Garlsbad, GA, USA). Sequencing was 
performed by a genetic analyzer (Apphed Biosystems 3130 
Genetic Analyzer, Garlsbad, GA, USA). Editing of the obtained 
sequences was carried out using the BioEdit software package 
(http://jwbrown.mbio.ncsu.edu/Bio-Edit/bioedit.html). Primer 
sequences (T3, T7, 571, 898E) were removed from sequences, 
then translated to protein sequences, and compared to known 
sequences using BLAST. Protein sequence alignments and 
phylogenetic analysis was done in MEGA V [34]. Sequences 
have been deposited in the GenBank database under accession 
numbers JQ6592 14 to JQ659253. 

Results 

In spring 2005, both sites (RK05 and ZK05) were studied in an 
initial test to determine if chemoautotrophy rates could be 
quantified by ''^G-DIG labeling of PLFA in the dark. Sites were 
sampled again in spring 2006 (RK06) and autumn 2007 (ZK07), 
when a more extensive sampling program was executed. 

Sediment biogeochemistry 

Oxygen penetrated significanfly deeper in ZK sediment (1— 
2 mm) than RK sediment (0.2-0.5 mm, Table 1). At RK06, high 
concentrations of free sulfide were found in the very top layer of 
the sediment, whereas sulfide only started to accumulate below 
2 cm sediment depth at ZK07 (Fig. 1). Sulfide and ammonium 
concentrations were more than 10 times higher for RK06 than for 
ZK07 throughout the sediment column (Fig. 1). In 2005, 
porewater samples were also taken at both sites and analyzed for 
suffide, but samples were taken two weeks before chemoautotro- 
phy measurements and not at exactiy the same location, which is 
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especially important for the RK site due to its patcliy nature. 
However, the contrast between the RK and ZK sites was similar 
with high concentrations of free sulfide at RK05 right to the top of 
the sediment core and no detectable sulfide in top 5 cm of the 
ZK05 sediment (results not shown). Porewater concentrations of 
Die and S04^' showed little variation with depth for ZK()7 
strongly indicating bio-irrigation, whereas DIG increased and 
S04^" decreased with depth for RK06 indicating carbon 
mineralization by sulfate reduction (Fig. 1). 

Diffusive sediment oxygen consumption rates as determined 
from microelectrode profiles were \'cry high for RK06 with 
192 mmol m~^ d~' and were approximately 15 mmol m~^ d~' 
for all other samplings (Table 1). The difference between the two 
RK samplings is probably due to the patchy nature of the site, 
even though visually similar black sediments with a whitish top 
layer were sampled in both years. For RK06, anaerobic carbon 
mineralization rates were about twice as high in the top centimeter 
(6.8 ±0.5 |J,mol C cm ^ d ') than in the 1—5 cm layer 
(3.2 ±0.8 (tmol C cm^'' d~'), whereas both sediment layers 
showed similar carbon mineralization rates for ZK07 (0-1 cm, 
1.6±0.9 nmol C cm"'^ d"'; 1-5 cm, 2.3±0.3 nmol C cm"^ d"'). 
Integrated over the upper 5 cm, anaerobic carbon mineralization 
rates were 197 and 106 mmol C m"^ d"' at RK06 and ZK07, 
respectively (Table 1). Ammonium production-based mineraliza- 
tion rates agreed with carbon mineralization rates at all sediment 
depths and ammonium-based mineralization rates in the top layer 
(0-1 cm) were 8.0±2.0 and 1.1±0.3 jtmol C cm"^ d"' for RK06 
and ZK07, respectively. 

Chemoautotrophy rates 

Dynamics of PLFA labeling with '^C-bicarbonate were studied 
in detail for RK06. Substantial labeling could already be detected 
in the 0-0.5 cm horizon after 4 hours of incubation and, although 
there was some variation, total PLFA labeling increased linearly 
with time for up to 4 days (R^ = 0.77, n = 8). Similar results were 
obtained for RK05 and ZK05 as calculated chemoautrophy rates 
were similar after 2 and 4 days of incubation (data not shown). For 
RK06. the "C- enrichment in the DIG pool in the 0—0.5 cm of the 
sediment decreased from 1800± 120 %o A8''^G (''^G enrichment of 
the DIG pool of 1.9"/o) approximately after 4 hours to 550 ± 1 10 %o 
AS'^^G (0.5% '"^G) after four days, probably because of exchange 
with atmospheric carbon dioxide and dilution with DIG produced 
during organic matter mineralization. Gores from RK06 could not 
be kept closed at the top because sub-oxic conditions developed 
within one day due to the very high oxygen consumption rates. 
The reported chemoautotrophy rates have been corrected for this 
change in DIG labeling with incubation time. 

Chemoautrophy was generally limited to the top centimeter of 
the sediment especially at the RK site (Fig. 2). For RK06, the main 
activity (95 ±1%) was found in the top 0-0.5 cm horizon, which 
contains the oxic top layer, and below 1 cm depth no activity 
could be detected. Similar results were obtained for RK05 when 
depth profiles were determined after 2 and 4 days. Interestingly, 
chemoautotrophy rates recorded in the top layer of RK06 were 
similar (7.2 Jtmol G cm ' d \ Fig 2) to the anaerobic carbon 
mineralization rates (6.8 )j.mol G cm~^ d~', see above), suggesting 
balanced GO2 production and consumption. For the ZK 
sediment, highest chemoautotrophy rates were also found in the 
top layer of the sediment, but activity remained relatively high 
down to 2 cm depth for both sampling dates. 

Depth integrated (0-5 cm), whole-core chemoautotrophy rates 
ranged from 3-36 mmol G m ^ d ' (Table 1). Rates measured at 
RK06 were very high, about 6 times higher than for RK05, which 
is probably due to the patchy nature of the site and differences 



between years (Table 1). Rates for ZK05 and ZK07 were however 
in the same range and always lower than for the RK site. Whole- 
core chemoautotrophy rates generally scaled with diffusive oxygen 
uptake rates and whole sediment chemoautotrophy to oxygen 
consumption ratios were relatively similar for all sites ranging from 
0.17 to 0.32 mol G (mol Oj)"' (Table 1). 

Active chemoautotrophic bacterial communities 

As we determined '''G-dark fixation rates into biomarker PLFA, 
the results can also be used to describe and compare active 
communities of chemoautotrophic bacteria. As an example, Fig. 3 
shows the PLFA concentration and labeling data as obtained 
during this study after one day of incubation for the top layer (0- 
0.5 cm) of RK06. The PLFA detected were typical for intertidal 
sediments with high amounts of both bacterial and eukaryote- 
specific biomarkers; the latter were probably mainly derived from 
the diatoms growing on the sediment surface as they were 
dominated by 20:4(n6 and 20:5(n3 that occur in high amounts in 
diatoms (Fig. 3A; [35]). Detection of PLFA labeling is based on the 
increase in 5''^G ratios and these were well above detection limits 
in many PLFA (Fig. 3B; 0.6 to 2%o A5'^C detection limit 
depending on compound). The AS'^'G ratios were highly variable 
between PLFA ranging for instance from 0 to 1 10%o after one day 
(Fig. 3B), suggesting that a specific sub-group of the total bacterial 
community was active. At RK06, PLFA that gained most ' G label 
were 14:0, 16:1<b7c, 16:l(n5, 16:0 and 18:1(b7c together explaining 
83 ±4% of the total incorporation into PLFA (Fig. 3G). There 
were also minor amounts of label recovered in 14:1, 15:1, 15:0, 
17:1(d8 and cyl7:0. Branched, i- and a-PLFA and eukaryote 
PLFA like 20:4<n6 and 20:5(03 gained very little label even though 
they were a dominant feature in the PLFA concentration pattern 
(Fig. 3A & 3G). 

The PLFA labeling pattern for the two ZK samplings was very 
different from the RK site strongly suggesting that difierent 
chemoautotrophic communities were active at the two sites. For 
comparison, the labeling pattern of ZK05 is also presented in 
Fig. 3D as it showed the largest differences from RK06. Main 
differences were a much higher labeling in all branched i- and a- 
PLFA and in several mono-unsaturated PLFA (15:1, 16:la)7t, 
17:lro8c, 18:lco9c and 18:l(n5) and uneven numbered saturated 
PLFA (15:0, {yl7:0 and 17:0) for both ZK samplings. 

To study the differences in PLFA labeling patterns further we 
performed a PGA analysis for all sediment layers where significant 
chemoautotrophy was detected (depth ranges RK 0-1 and ZK 0- 
2 cm; Fig. 4). The first PGA axis explained 42% of the variation 
found in the data set, whereas the second axis added another 1 5 % . 
Glustering was mainly based on I6:1(b7c versus 15:0, 17:0, 
17:l(n8c, 18:0 and all branched PLFA for tlu- first axis similarly as 
seen in Fig. 3 and 16:0 and 16:lco7t versus 14:0, 14:1 and cyl7:0 
for the second axis. Glustering of sediment samples was mosdy 
determined by sampling site and sampling year with both RK 
samplings clustering closely together and more dispersal amongst 
the samples from the ZK site was observed (Fig. 4A). Some 
additional variation was found with sediment depth but only for 
ZK06, where the top layer data (0-1 cm) clustered more closely 
together with the RK data whereas the deeper layers were shifted 
towards the ZK05 samples (Fig. 4A). Distribution of label among 
PLFA did not change substantially with incubation time, 
suggesting that active chemoautotrophic communities remained 
similar for up to 4 days. The PGA analysis therefore also showed 
major difference in active chemoautotrophic bacterial communi- 
ties between the high free-sulfide RK and low free-sulfide ZK 
sediments. 
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Figure 1 . Porewater concentrations for the RK06 and ZK07 sediments. Shown are data for sulfide and ammonium (A, B) and for sulfate and 
Die (C, D). 

doi:1 0.1 371/journal.pone.01 01 443.g001 



There was also indirect evidence of transfer of dark-fixed carbon 
to fauna: some ' ^C-labelling was detected in bulk sediment PLFA 
characteristic of Eukarya and therefore fauna, 18:2w6c for both 
sites and 20:4(b6 and 20:5co3 for die ZK site (Fig 3C and D). 

Rubisco type lA diversity 

To further characterize the chemoautotrophic community, we 
used a novel primer set to construct clone hbraries for the Rubisco 
Type lA large-subunit gene for both sites in 2008. For the RX site, 
1 7 clones related to Type lA Rubisco were recovered and site ZK 
yielded 23 clones. We also found a limited number of clones 
related to unicellular cyanobacteria derived type IB Rubisco 
especially for the RK site (not shown). 

The phylogenetic relationship for the clones from l)()th sites 
together with other environmental clones and related micro- 
organisms is shown in Fig. 5. Type lA clones were found in two 
main clades labeled type lA 1 and lA 2 (Fig. 5). Rubisco type lA 1 
clones were most closely related to uncultured faunal endosym- 
bionts belonging to the Gammaproteobacteria and marine 



sediment clones from a variety of other studies. Whereas sequences 
in clade lA 2 were- related to various chemoautotrophic sulfur and 
ammonium oxidizing bacteria and also to other environmental 
sediment clones. Be^iatoa Rubisco type LA also clustered in clade 
lA 2, but although clones from both RK and ZK were found in 
this clade they were not closely related to Beggiatoa suggesting that 
they belonged to other groups of chemoautotrophic bacteria. In 
general rather similar sequences were recovered from both sites, 
although ZK clones were relatively more abundant in clade LA 1 
and RK clones dominated clade LA 2. The restdts suggest that 
chemoautotrophic communities, which use the Calvin cycle for 
carbon fixation, were relatively similar at both sites and that 
Beggiatoa could not be detected at least for the 2008 sampling. 

Discussion 

Rates of chemoautotrophic carbon fixation 

We detected very high rates of chemoautotrophic dark fixation 
in the top layers of two intertidal sediments especially for the 
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RK06 sediment. Volumetric chemoautotrophy rates detected in 
the top of the RK sediment were comparable to some of the 
highest sulfate reduction rates detected in marine sediments [36]. 
These high rat(-s arc in part explained because reduced substrates 
produced by mineralization processes are released throughout the 
active sediment column whereas the reoxidation by chemoauto- 
trophic bacteria is concentrated in the top of the sediment. Dark 
fixation rates as measured in our study also includes anapleurotic 
reactions by heterotrophic bacteria, which may account for up to 5 
to 10% of the biomass produced by in all bacteria including 
heterotrophs [37-39]. However, chemoautotrophy rates detected 
in the top-layer of the RK06 sediment were actually similar to 
carbon mineralization rates (Fig. 2 and see Results), suggesting 
that the role of anapleurotic reactions was minimal explaining at 
most about 5 % of the measured dark-frxation rates if one assumes 
a relatively high heterotrophic growth efficiency of 50% [40]. An 
additional advantage of measuring dark fixation through '^C- 
labehng of PLFA is that carbon fixed through anaplerotic 
reactions is not directiy utilized in the synthesis of lipids such as 
fatty acids [39,41]. Based on oxygen consumption rates and 
assuming a coupled system, chemoautotrophy explained between 
18% (RK06) and 32% (RK05) of die sediment carbon cycling, 
which would make chemoautotrophy the second or third most 
important biological carbon cycling process after anaerobic carbon 
mineralization and possibly photosynthesis by benthic diatoms in 
these intertidal sediments. 



Sulfur oxidizing bacteria are expected to be the main 
chemoautotrophs in these intertidal sediments; the contribution 
from ammonium oxidizing bacteria should be less important 
because about 6 timc-s more sulfide than ammonium is produced 
during anaerobic carbon mineralization given the typical C:N 
ratio for marine organic matter [42]. In addition, nitrifying 
prokaryotes also tend to have lower growth ^delds per mol of 
substrate oxidized than sulfur oxidizers [43]. The chemoautotro- 
phy data from our study scaled well with measured diffusive 
oxygen fluxes with an average whole system yield (± SD) of 
0.22 ±0.07 mol C (mol Oj)"' (Table 1), which is very similar to 
the typically reported maximum growth yields for aerobic sulfur 
oxidizing bacteria of 0.23±0.11 mol C (mol Og)"' [44-50]. The 
similarity in C:02 yield between intc'rtidal sediments and sulfur 
oxidizing bacterial cultures can only be explained if most of the 
sediment oxygen cxmsumption was indeed used for reoxidation of 
reduced sulfur and if reoxidation was predominantly performed by 
obligate sulfur-oxidizing chemoautotrophic bacteria growing close 
to their maximum reported yields. In addition, chemoautotrophic 
bacteria should have efiectively competed with chemical oxidation 
processes such as the oxidation of free sulfide with oxygen, as has 
been shown in gradient systems where oxygen and sulfide are 
found in close proximity similar to the RK site [3]. Furthermore, 
even though these coastal sediments receive very high organic 
matter inputs, our data suggest that the activity by heterotrophic 
and mixotrophic sulfur-oxidizing bacteria was also limited as this 
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should have led to lower C:02 yields. Competition between 
autotrophic and heterotrophic or mixotrophic sulfur oxidizing 
bacteria depends to a large extent on the sulfur/ organic substrate 
ratio with low ratios supporting heterotrophic sulfur oxidation 
[11,51], but this apparendy didn't play a role in the studied 
sediments possibly due to strong competition for organic substrates 
by other heterotrophic bacteria. Our re.sults suggest that dark 
fixation rates as determined by ' ^C-bicarbonate labeling of PLFA 



yield realistic chemoautotrophy rates in relation to sediment 
biogeochemistry and that chemoautotrophic bacteria were grow- 
ing with high efficiency independent of sediment sulfur chemistry. 

Chemoautotrophy rate data are available from only four other 
studies on coastal marine sediments, all of which measured total 
dark fixation rates by determining '*C-DIC incorporation into 
POC [12-15]. It should be noted that two of these studies are 
based on laboratory incubations with homogenized sediments 
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[12,15]. Thomsen and Kristensen [12] reported maximum rates of 
approximately 0.35 |imol C cm~'^ d~' (averaged over the same 
depth as in the present study) similar to site ZK, whereas rates 
found by Enoksson and Samuelsson [13] and Lank at al [14] are 
lower at about 0.12 (imol C cm ^ d '. The data collected by 
Thomsen and Kristensen [12] yield a C:02 ratio of 0.24 mol C 
(mol O.)) ' similar to our data, whereas the chemoautrophy rates 
reported by Enoksson and Samuelsson [13] are substantially 
higher than expected from their oxygen consumption rates (C:02 
ratio about 1.2 mol C (mol 02) '). The C:02 ratios from 
Enoksson and Samuelsson [13] can not be readily explained in 
relation to the reported growth yields of sulfur reducing bacteria 
[44-50], suggesting that the chemoautotrophy rates for this 
sediment may have been substantially overestimated possibly due 
to the incomplete removal of '''C-DIC label. Santoro et al [15] 
found much lower yields of about 0.025 mol C (mol 02) ' for 
three brackish coastal lake sediments. Lenk et al [14] also reported 
C-POC based chemoautotrophy rates of 4. 16±0.03 mmol C 
m ^ d ' for an intertidal flat consisting of permeable sands, which 
related well to sediment sulfide fluxes [14,52]. Lenk et al [14] did 
not report oxygen consumption rates. However, potential oxygen 
fluxes of approximately 70 mmol O2 mT ^ d~' have been reported 
for the same sand flat [53], suggesting a relatively low C:02 yield 
of about 0.06 mol C (mol 02)"'. Oxygen fluxes may have been 
overestimated because they are potential rates based on aerobic 
incubations with sediments from different horizons that may not 
always be oxic. The depth distributions of chemoautotrophy of 
Enoksson and Samuelsson [13] and Lenk et al. [14] were diflerent 
from our study as they both found substantial rates deeper in the 
sediment up to a depth of 10 cm. Lenk et al [14] studied a 
permeable sediment where oxidants such as oxygen are trans- 
ported deep into the sediment by advective porewater flows [53], 
which may explain the high chemoautotrophy rates deeper in the 
sediment. The chemoautotrophy rates in our study at the RK06 
(3.2 to 6.8 \ixao\ C cm"'* d~') are the highest reported so far for 
coastal marine sediments. 



Furthermore, our results suggest that chemoautotropliically 
fixed carbon may also be an important food source in the 
microbial food web in these typical coastal sediments. Chemoau- 
totrophy rates were similar to carbon mineralization rates in the 
top layer of the RX06 sediment (Fig. 2 and Results). Net 
consumption of DIG related to the chemoautotrophy has been 
indicated in the top layer of active coastal sediments [12,54]. This 
suggests that the production by chemoautotrophs dominated the 
microbial food web and that heterotrophic bacterial secondary 
production was less important (if one assumes a growth efficiency 
of 50% for heterotrophic bacteria). Santoro et al [15] measured 
growth of both chemoautotrophic and heterotrophic bacteria in 
three brackish lakes and found that chemoautotrophy could 
explain up to 50% of the heterotrophic bacterial growth. 
Additionally, chemosyntheticaUy produced biomass may poten- 
tially be an important source of energy fueling the benthic food 
web. We indeed detected limited labeling in eukaryotic fauna- 
derived PLFA for both sediments (Fig 3C & D), which suggests 
that sediment fauna may in part be feeding on chemoautotrophic 
bacteria. Chemoautotrophic bacteria support the food web in 
many extreme marine ecosystems like hydrothermal vents and 
mud volcanoes with limited organic matter inputs [8,55] . Based on 
our results, the role of chemoautotrophic carbon fixation in the 
benthic food web of coastal marine sediments should receive more 
attention. 

Communities of chemoautotrophic bacteria 

We detected clear difference in PLFA labeling patterns for the 
ZK and RK sediments suggesting that the active chemoautotro- 
phic bacterial communities were substantially different (Fig. 3 and 
4). The classical sulfur and ammonium oxidizing bacteria 
predominantly contain even-numbered saturated and mono- 
unsaturated PLFA such as 14:0, 16:1(b7c, 16:0 and 18:l(jo7c 
[17,55-60]. For the RK sediment, most of the label was indeed 
recovered in these typical PLFA and the PLFA labeling pattern for 
instance closely resembles the PLFA composition of Beg^atoa 
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Figure 5. Phylogenetic relationships between Rubisco Type 1 A clones recovered from the high f ree-sulfide RK and low f ree-sulfide 
ZK site in 2008 together with relevant sequences of environmental clones from other studies and chemoautotrophic prokaryotes 
retrieved from Genbanit (El: [73]; NK: [32]; Gi: [75]; Md (Madrid et al. unpublished)). The Neighbor-Joining tree is based on amino acid 
sequence and bootstrap values are based on 1000 times replication. 
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especially for RK()5 (Fig. 3 and [55,56]). Some types of sulfur 
oxidizing thiohacilli also contain uneven numbered PLFA like 
cyl7:0 [17,57], and label was indeed recovered in these PLFA 
especially at both ZK samplings and at RK06 suggesting that these 
groups may also be active in coastal sediments. In view of the high 
free sulfide and ammonium cxmcentrations in the porewater, the 
PLFA labeling at the RK site was therefore in agreement with the 
activity of typical chemoautotrophic sulfur and ammonium 
oxidizing bacteria. 

In contrast, major amounts of label were also found in branched 
i- and a-PLFA at the ZK site, which have not been reported in the 
classical chemoautotrophic sulfur oxidizing and nitrifying bacteria. 
However, several groups of sulfate reducing bacteria contain large 
amounts of branched PLFA [61-63]. In addition, significant 
amounts of label were recovered in il7:l(n7, al7:lco7, 17:1cd8c 
and cyl7:0/17:lcB6c that have been suggested as specific 
biomarkers for certain groups of sulfate reducing bacteria [61- 
63], suggesting that sulfate reducing bacteria may be important 
chemoautotrophs in the ZK sediment. Chemoautotrophic growth 
can be found in two groups of sulfate reducing bacteria namely 
those that use hydrogen gas as substrate [64] and sulfur 
disproportionating bacteria that gain energy from inorganic 
fermentation of substrates such as thiosulfate and elemental sulfur 
to sulfate and sulfide [65,66]. Both hydrogen turnover and 
disproportionation reactions have been shown to be important 
processes in anaerobic marine sediments [67,68] and they have 
been indicated to possibly support chemoautotrophy in coastal 
sediments [12]. For the ZK07 sediment, the labeling in these 
branched PLFA became progressively more important with depth 
(Fig. 4), suggesting that it may indeed be associated with these 
anaerobic processes. Bran[:hed PLFA have also been found in 
anammox bacteria [69], but these also contain 10Mel6:0 which 
showed no or Umited labeling in our study (Fig. 3) suggesting that 
they were not important as chemoautotrophs. This is in agreement 
with studies that found that denitrification and not anammox is the 
dominant nitrate consuming process in active coastal sediments 
[70,71]. The high labeling in branched PLFA as detected at the 
low free-sulfide ZK site therefore indicates that anaerobic 
chemoautotrophy most likely by bacteria related to sulfate 
reducers may be important in coastal sediments. 

We also studied the diversity of chemoautotrophic bacteria that 
use the Calvin cycle for carbon fixation by constructing clone 
Ubraries of the Rubisco type lA gene. Clones detected fell in two 
main clades: one related to uncultured endosymbiontic sulfur- 
oxidizing Gammaproteobacteria (clade lA 1) and another related 
to free living chemoautotrophic sulfur and ammonium oxidizing 
bacteria (clade lA 2, Fig. 5). Although clade lA 1 sequences were 
most closely related to endosymbionts, they are also commonly 
found in related uncultured free-living i)acteria [14,72]. These two 
Rubisco clades are usually also detected in other environmental 
studies on marine sediments [32,73]. Lenk et al [14] also showed 
that that free-living Gammaproteobacteria related to sulfur- 
oxidizing endosymbionts were important for dark carbon fixation 
in a permeable tidal flat sediment. The differences in Rubisco 
clone distribution between the two sediments are most likely due to 
differences in the availability of reduced sulfur compounds namely 



free sulfide at the RK site and iron sulfides at the ZK site. 
However, it seems unlikely that these differences in clone 
distribution explain the contrasting PLFA labeling patterns from 
the two sites. First, the differences in Rubisco clone distribution 
appear to be less site-dependent compared to the PLFA labeling 
patterns that showed major differences between sites. Second, the 
detected clones are related to bacteria that typically do not contain 
i- and a-branched PLFA as found for the chemoautotrophs from 
the ZK site. As discussed before, the PLFA composition of both 
the typical chemoautotrophic sulfur and ammonium oxidizers in 
clade lA 2 and Gammaproteobacteria like the ones detected in 
clade lA 1 is usually dominated by straight saturated and un- 
saturated compounds similar to the labeling patterns at the RK 
site [17,57,58]. Third, sulfate reducing bacteria typically fix 
inorganic carbon by either the reversed TCA-cycle or the 
reductive acetyl-CoA pathway [64] and are therefore not targeted 
in our Rubisco assay. The dominance of Gammaproteobacteria 
among the Rubisco clones is therefore also agreement with the 
PLFA labeling patterns especially for RK and to lesser degree for 
ZK. 

To conclude, we detected ver\' high rates of chemoautotrophy 
in two active intertidal sediments and our results indicate that 
chemoautotrophic carbon fixation is an important part of the 
carbon cycle of coastal sediments. Clear differences were found in 
active chemoautotrophic bacterial communities probably caused 
by differences in reduced sulfur compounds available in the two 
sediments. Our study suggests that anaerobic bacteria related to 
sulfate reducers played an important role at the low-sulfide ZK site 
while typical sulfur-oxidizers, probably Gammaproteobacteria, 
were more important at the free-sulfide RK site. Chemoautotro- 
phic Archaea such as the ammonium-oxidizing Crenarcheoata 
widely found in marine sediments [74] were not considered in the 
present study, but their activity could be studied by determining 
'^C labeling in Archaeal ether-lipids [41]. The number of coastal 
sediments where chemoautotrophic carbon fixation rates have 
been determined is still very low especially compared to studies on 
other aspects of the sedimentary carbon cycle. Besides our study 
we are aware of only four other studies where chemoautotrophy 
rates were quantified in coastal marine and brackish lake 
sediments [12-15]. Hence, chemoautotrophy in coastal marine 
sediments as driven by reoxidation processes should receive more 
attention in future studies given its importance in the carbon cycle 
of coastal sediments and the potential role in the benthic food web. 
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